Theory of Neutron Diffraction from the Vortex Lattice in UPt3 
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Abstract 

a 

a. 

. Neutron scattering experiments have recently been performed in the super- 

C 

conducting state of UPt3 to determine the structure of the vortex lattice. 

O 

The data show anomalous field dependence of the aspect ratio of the unit cell 
in the B phase. There is apparently also a change in the effective coherence 
length on the transition from the B to the C phases. Such observations are not 
consistent with conventional super conduct vity. A theory of these results is 
constructed based on a picture of two-component superconductivity for UPt3. 
In this way these unusual observations can be understood. There is a possible 
discrepancy between theory and experiment in the detailed field dependence 
of the aspect ratio. 
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Experiments of the last several years on superconducting UPt3 have revealed a rich phase 
diagram. The phase boundaries have been mapped with considerable accuracy by locating 
anomalies in specific heat []J [Q, ultrasound ||, and torsional oscillator signals The 
additional transitions found may be explained by assuming that the superconductivity is 
unconventional, and that the order parameter has two complex components. At zero 
pressure, there are three superconducting phases, called the A phase (low fields and high 
temperatures), the B phase (low field, low temperature), and the C phase (high field, low 
temperature). The theoretical explanation remains untested, however, unless detailed ex- 
perimental characterization of the phases can be carried out, preferably by searching for 
definite signatures of the predicted phases. The theory suggests that the B phase is the 
most promising in this regard. In this phase, the two components are predicted to coexist, 
and their interplay can give rise to novel effects. 

The comparison of theory and experiment is not straightforward, however, as no micro- 
scopic experimental probe couples directly to the order parameter. Because of this, it is 
desirable to have a probe for which predictions can be made using the phenomenological 
Ginzburg-Landau theory alone, thus circumventing microscopic details. Neutron diffraction 
from the flux lattice depends only on the field distribution in the sample |J and therefore 
may be calculated from the phenomenological theory, yet it offers a level of detail which 
cannot be obtained by any other means 0. A new experiment of this kind has recently 
been carried out and is reported in the accompanying Letter 0. 

Calculations of the lattice structure and the field distribution have so far been limited 
to the neighborhood of the upper critical field H c2 (K| However, it is only at low 



temperatures that the penetration depth is short enough that experiments can be done. 
In practice this means well away from H c2 . Furthermore, the interesting effects occur in 
the field dependence of the scattering intensities in the B phase, and particularly near H*, 
the transition field from the B to the C phase. The purpose of this Letter is to present 
calculations of the lattice structure and the scattering cross section as a function of applied 
field, and to compare the results to experiment. 



The scattering amplitude for a change in neutron wavevector Q in the Born approxima- 
tion is H 

MQ) = ' mN - N I d 2 rh(r)e^ = h(Q), (1) 

where and fijy are the mass and magnetic moment of the neutron and ki is the incoming 
wavevector. h(r) is the field in the sample and h(Q) is its Fourier transform. The integration 
is over a cross-section A of the sample perpendicular to the beam. The object of interest, 
the observed intensity I{Q), is thus proportional to \h(Q)\ 2 . The experiments have all 
been performed with H, the applied field, directed in the basal plane, so the calculations 
presented here will be only for that field direction. I deal only with Bragg scattering, so Q 
is a reciprocal lattice vector. 

If the y-axis (in the basal plane of the hexagonal crystal) is taken as the field direction, 
two-component theories have the generic free energy density: 

f(ff) = a(\ Vx \ 2 + \ Vy \ 2 ) + fitf- rf) 2 + m- V\ 2 (2) 
+K x \D xVx \ 2 + K y \D xVy \ 2 + K z (\D zVx \ 2 + \D zVy \ 2 ). (3) 

7] x and rjy are the two components of the order parameter. D x = — id x — 2eA x /hc and 
similarly for D z . r\ x and r\ y may belong to one of the two-dimensional representations E\ or 
E 2 of the point group, or they may belong to different, but accidentally nearly degenerate 
representations. The scattering experiments probably do not distinguish these possibilities. 
For simplicity, I have taken t] x and f] y to be degenerate (only a single a) at zero field. This 
is reasonable at the low temperature (50 mK) at which the experiments are done. In fact I 
shall fix the temperature and treat a as constant. This means that only deal with the B (low 
field or H < H*) and C (high field or H > H*) phases, and the transition between them, 
will be discussed. The free energy F = J f[fj(r)]d 3 r where f[ff(r)} is given by Eq. |3| does 
indeed lead to the observed phase diagram with three superconducting phases for H in the 
basal plane [[J. The form of rj(r) in the C phase is known. It is found that rj(f) = (i] x (f), 0). 
Thus the free energy reduces to 
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fc(v) = a\ Vx \ 2 + (A + (3 2 )\ Vx \ A + K x \D xVx \ 2 + K z \D zVx \ 2 . (4) 
The problem represented by this free energy density is isomorphic to that of an s-wave 



superconductor with mass anisotropy. One may therefore transcribe well-known results |L2 
|*tf>|. = * [1 + + (Af) 2 QS]- 2 exp {-\ [ ( &Ql + «? )»«»]} , (5, 

where £p = (JQ/a:) 1 / 2 are the coherence lengths and Af = (7i 2 c 2 /327re 2 lfi|?7 a ;| 2 ) 1 / 2 are the 
penetration depths in the C phase. Since the free energy density may be transformed 
to the isotropic form by rescaling coordinates: x' = (£, z /£x)^ 2 x and z' = (^/^) 1//2 ^, the 
fluxons form a rescaled hexagonal lattice (centered rectangular lattice). In London theory the 
vortex cores act as delta-function sources of the field. Corrections to this require numerical 



calculations which have been carried out by Brandt ||13|| , and it is concluded that the sources 
are well represented by Gaussians - this is reason for the exponential factor in Eq. [|. The 
coordinate scaling determines the opening angle ai of the reciprocal lattice (defined in Fig. 
1), which uniquely determines the aspect ratio of the unit cell. It is given by tan 2 (a^) = 
3K X / K z . This result and flux quantization allow us to calculate the reciprocal lattice vectors 
as a function of field. Substituting these values into Eq. | shows that I(Q) is the same for 
all Q in the first shell (the six smallest nonzero Q). These are the only points measured to 
date. 

Theory thus predicts that three properties qualitatively characterize the C phase: 

(i) the lattice structure (shape of the unit cell) is independent of field; 

(ii) all peaks in the first shell have the same intensity; 

(iii) the intensities fall off exponentially with field, with —dlnl/dH ~ ■ 

Property (i) has been pointed out before fl4jl , [fTD| . The C phase has no unique signatures 
of unconventional superconductivity, however, since an s-wave superconductor with mass 
anisotropy has all of these features. 

The B phase is quite different. I first develop the theory and then turn to comparison 
with experiment. 
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The opening angle in the B phase can be computed in the regime where £j << a « Xi, 
where a is the lattice constant for the vortices. In terms of the field, this is H cl << H << 
H C 2. Since UPt3 is strongly type-II, this is a substantial range. It will suffice for our purposes 
to compute the currents at a distance of order a from the cores, because the momentum 
transfers of experimental interest are of order 1/a. Thus the structure of the cores at short 
distances of order a complex subject, is not of interest here. At these larger distances 
the /?2 term in Eq. || locks the relative phase of r\ x and r\ y : r\ y = ±ir(H)i] x , where r(H) is a 
real ratio. r(H) supplies the interesting field dependence in the B phase. Since ff y appears 
continuously, r(H) is a nonnegative, monotonically decreasing function of H with r(H*) = 0. 
The phase-locking relation, together with the free energy of Eq. [3], gives a London equation 
for the currents with the penetration depths 

h 2 r 2 

(\B\2 = '±± ( 6) 

1 x) 32ire 2 \r) x \ 2 [K x + r 2 (H)K y Y 1 ; 

and 

h 2 c 2 

(\B\2 _ 'J_± ( 7 \ 

V z) 32ire 2 \7] x \ 2 K z [l + r 2 {H)Y 1 ' 

In the field regime under consideration, these currents determine the lattice structure and 
one finds 

2 K x + r 2 (H)K y 

tan K) = 3[1 + r 2 (H)]K z - (8) 

The shape of the unit cell is strongly field- dependent in the B phase. This effect does not 
occur in one-component superconductors, conventional or unconventional. 

The computation of the intensities in the B phase is not so straightforward, since cor- 
rections to London theory are involved. However, the same phase-locking approximation 
allows us to make a mean-field-type theory. We neglect amplitude correlations in the region 
where the distance from the cores is much greater than the coherence lengths. Then the 
effective field on either component has the same spatial dependence as in the s-wave case, 
and we may again apply the results of Brandt. The complication which arises is that both 



rf x and fj y act as sources of the field. This leads to separate exponential dependences, and 
the breakdown of the simple relationship —dlnl/dH ~ £ 2 . There are now two effects which 
determine the field dependence of the intensity. First, the interaction of ff x and ff y given by 
the quartic terms in Eq. |3] sets up effective fields which give a field dependent renormaliza- 
tion of the correlation lengths of both components. This means that the exponents acquire 
an additional field dependence relative to Eq. |[ Second, the prefactors of the exponents 
have a field dependence because of the separate contributions from ff x and ff y , whose relative 
weight is field-dependent. 

The full expression for the h(Q) may be separated into a part which depends relatively 
weakly on H and is probably unobservable, and the exponential part The full expression, 
and details of its derivation, will be given elsewhere. The exponential part is 



/'(t?V-[(i + /- 2 r J F(0,^Q. i eQ,) (9) 
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+(i + t'kjk^f^Cq^CM] x ex V [(-(C?Ql - (£?Ql)/q (io) 
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+ (r- 2 + KjK y )-'F{^t%Q z ,i* y Q x )] x exp[(-(^) 2 Q^ - (12) 
In this formula, the coherence length t^ x is given by 

& - <f > i/2 1 1 + > - 1»< 2 - h - 1» - M' 2 • (i3) 



and there are similar formulas for the other three coherence lengths. The field H y is a 
constant given by 

H* 

Hy = l-b + bH'/Ha' (14) 
with b={(3\ — ^2) /{Pi + ^2); while the function F(x, gi, q 2 ) is defined as 

/ — oo r — 00 ^2 

du / dv — -exp[-(u-iq 1 /2) 2 - (v -iq 2 /2) 2 }. (15) 
-CO J — CO X U ~r V 
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Thus the effective coherence length, if it is defined by the slope of In I(Q) as a function of 
field, is seen to be field dependent. This again is utterly characteristic of a multicomponent 
system, and cannot be found in a one-component superconductor. 

These calculations have been done assuming that the separation of the singularities in 
ij x and fj y in the unit cell are separated by a distance much less than the lattice constant. If 
this is not the case, as in the 'shifted' phase predicted near the tetracritical point, then the 
intensities for the reciprocal lattice vectors in the first shell may be quite different because 
of extinction effects. In particular, at low field, the intensity for Q ,i will be less than that 
for Qi-\. Unfortunately, explicit calculations for this phase are quite difficult. 



Eqs. |8| and [12] give three qualitative predictions for the B phase: 

(i) the lattice structure (shape of the unit cell) is depends on field; 

(ii) peaks in the first shell have different intensities; 

(iii) the intensities have a complex field dependence, with deviations from the exponential 
shape; 

(iv) there is a kink in 1(H) at the second-order B-C transition. Property (i) has been pointed 
out before 0, fH] , fllOfl . Property (iv) does not imply that the BC transition is first-order. 
The slope of 1(H) is not interpreted as proportional to a single coherence length, (which 
would then be discontinuous). Instead, the kink is interpreted as signalling the continuous 
growth of a new component of the scattering. The BC transition is second-order in this 
theory. 

Before comparing these predictions quantitatively with experiment, we must discuss the 
determination of the Ginzburg-Landau parameters. The final results for the opening angle 
of the flux lattice depend only on the ratio P1/P2 an d the ratios of the stiffness constants 
K x , K y and K z . I take Pi/ (3 2 = 0.5, as determined by specific heat experiments. The stiffness 
constant ratios as determined by fitting to the neutron data are K x : K y : K z = 1.5 : 0.88 : 
2.5. 

The opening angle is plotted in Fig. 1. The theory predicts a field-dependent lattice 
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structure in the B phase: H < H* = 5.3kOe, as is observed. There should be a kink in 
the curve at H* and in the high-field C phase, should be independent of field. These 
predictions are consistent with the data, but more points at higher fields and smaller error 
bars are needed to provide a real test. 

The logarithm of the form factor Hi Q at the Bragg point Q\ t \ is plotted as a function of 
field in Fig. 2 with the same parameter ratios K x : K y : K z = 1.5 : 0.88 : 2.5. The coherence 
length, defined as the geometric mean of the coherence lengths in the x and z directions, 
was taken to be lOlA. Very good agreement with experiment is obtained. In particular, 
the truly novel feature in the data, the kink at H = H* = 5.3kOe, is very well reproduced. 
I have also calculated the intensity at the Bragg point Qo,i- It is not shown because for 
the parameter range here the calculated intensities differ by only a few per cent. It will be 
difficult to verify that this difference exists at the current level of experimental accuracy. 

The stiffness constant ratios obtained by the present fit can be compared to ratios ob- 
tained by fitting to the measured values of the upper critical fields at temperatures near 
T c for different field directions and from the discontinuity in slope at the tetracritical 
point. This gives K x : K y = 2, and applying constraints from particle-hole symmetry 
yields K x : K y : K z = 1:2:4. It is quite clear that the two methods disagree, and even the 
ordering differs between K x and K y . Indeed, if K x < K y , as suggested by the critical field 
slopes, the curve ai{H) is monotonically decreasing for H < H*, in clear contradiction to the 
data of Fig. (1). The significance of the discrepancy is not clear at present, since it involves 
an extrapolation from high to low temperatures which may not be justifiable. Nonlocal 
corrections at low temperatures are surely important and the nature of these corrections for 
unconventional superconductors is not known. 

The qualitatively new features predicted for neutron scattering from the flux lattice in 
the B phase are observed, as is the relatively conventional behavior of the C phase. Both 
the opening-angle data and the field-dependent intensities can be fit to good accuracy. This 
fit does not agree with one obtained from critical field measurements. This may be due to 
difficulties of extrapolation to the low temperatures of the experiment, or may indicate a 
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real discrepancy between theory and experiment. 
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FIGURES 

FIG. 1. Theoretical curve and experimental data points from Ref. 8 for the opening angle of 
the flux lattice as a function of the applied field. The inset shows the definition of the angle. 

FIG. 2. Theoretical curve and experimental data points from Ref. 8 for the logarithm of the 
form factor at the first Bragg peak as a function of the applied field. 
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